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Abstract 



We study carrier interaction induced many-body effects on the exci- 
tonic optical properties of highly photoexcited one-dimensional semiconduc- 
tor quantum wire systems by solving the dynamically screened Bethe-Salpeter 
equation using realistic Coulomb interaction between carriers. Including dy- 
namical screening effects in the electron/hole self-energy and in the electron- 
hole interaction vertex function, we find that the excitonic absorption is essen- 
tially peaked at a constant energy for a large range of photoexcitation density 
(n = — 6 x 10 5 cm -1 ), above which the absorption peak disappears without 
appreciable gain i.e., no exciton to free electron-hole plasma Mott transition 
is observed, in contrast to previous theoretical results but in agreement with 
recent experimental findings. This absence of gain (or the non-existence of a 
Mott transition) arises from the strong inelastic scattering by one-dimensional 
plasmons or charge density excitations, closely related to the non- Fermi liq- 
uid nature of one-dimensional systems. Our theoretical work demonstrates a 
transition or a crossover in one-dimensional photoexcited electron-hole system 
from an effective Fermi liquid behavior associated with a dilute gas of nonin- 
teracting excitons in the low density region (n < 10 5 cm -1 ) to a non- Fermi 
liquid in the high density region (n > 10 5 cm" 1 ). The conventional quasi- 
static approximation for this problem is also carried out to compare with 
the full dynamical results. Numerical results for exciton binding energy and 
absorption spectra are given as functions of carrier density and temperature. 

PACS numbers: 78.55.-m; 71.35.Cc; 78.66.Fd; 73.20.Dx; 71.10.P; 71.10.H 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



Excitons in low dimensional semiconductor systems have been extensively studied in 
the recent past. Present interest has focused on one-dimensional excitons in artificially 
structured semiconductor quantum wire (QWR) systems where spectacular improvements 
in growth and nanofabrication techniques have led to very narrow wires of nanostructure size 



conduction band of such a QWR most likely form a pure one-dimensional (ID) system. For 
the holes in the QWR valence band, the bare confinement potential (for example, in GaAs- 
AlGaAs system) is known to be too shallow (~ 10 meV) for a hole to be one-dimensionally 
confined in these QWR structures. Including the Coulomb interaction between electrons and 
holes along the transverse (i.e. perpendicular to the ID free motion direction) directions of 
the wire, however, Glutsch et al. pi] find that even the holes in the valence band of QWR 
can be strongly localized in the transverse plane, leading to both electrons and holes being 
effectively ID (or rather quasi-lD) in the dynamical sense. Therefore an exciton in such 
ultranarrow QWR nanostructures can be effectively thought of as a bound pair of a ID 
electron and a ID hole with the carrier dynamics being free along the ID wire direction as 
long as one is interested in low energy (lower than confinement energy ~ 20 — 100 meV) 
excitonic optical properties. Such strong confinement for both electrons and holes also 
substantially enhances the excitonic binding energy leading to novel optical phenomena. 
In the low (electron-hole) density limit without considering the self-energy correction to 
the conduction and the valence band energies as well as neglecting all dynamical screening 
effects, the single electron and single hole problem in forming the exciton can be exactly 
solved as a quasi- ID hydrogenic (Wannier exciton) atom with an exciton radius of about 

100 A for GaAs based QWR systems. This single ID exciton problem, where an electron 
and a hole in a QWR form a bound excitonic state, has been studied extensively in the 
recent literature in the context of understanding QWR excitonic optical properties. Such a 
non-interacting exciton picture, based on a simple single particle electron-hole hydrogenic 
bound state scenario, obviously only applies in the dilute low carrier density limit when the 
excitons or the bound electron-hole pairs are effectively very far from each other forming a 
noninteracting exciton gas. We will refer to this situation as a Fermi liquid (because in ID 
only an effectively noninteracting system may behave as a Fermi liquid) or a noninteracting 
exciton gas. In the high carrier density situation the excitons must overlap a great deal, 
and the simple Fermi liquid picture of a noninteracting exciton gas would not apply. Our 
main goal in this paper is to theoretically study this transition between the low density 
(Fermi liquid like) exciton gas and the high density system of interacting (and strongly 
overlapping) excitons in quasi-lD semiconductor (GaAs) QWR systems. This exciton gas 
to a strongly overlapping and highly correlated electron-hole system crossover with increasing 
electron-hole density can be thought of as a transition from an insulating exciton gas to a 
conducting electron- hole plasma (EHP), the Mott transition. A typical feature of this Mott 
transition, observed in higher dimensional (2D, 3D) optical experiments, is the development 
of optical gain in the absorption spectra where the absorption coefficient becomes negative 
(gain region) in some frequency range. One of the questions addressed in this paper is 
whether such an optical gain region exists in ID photoexcited QWR systems. In this paper 
we consider the formation, stability, and optical properties of one-dimensional excitons from 
low to high carrier densities in semiconductor QWR under photoexcitation conditions (i.e. 
equal electron and hole densities), a problem which has attracted a great deal of both 
theoretical [f^TH] and experimental attention in these years. Consistent with recent 
interest one of the central issues we focus on is the density-induced exciton gas to EHP Mott 
transition in ID QWR systems and its experimental signature. 

The motivation of our work arises from recent experimental studies of the photolumines- 
cence spectra of ID GaAs/Ali-^Ga^As semiconductor QWR systems jl|-fl]- The experimen- 
talists use strong lasers to pump photons into the QWR systems, exciting electrons from the 
filled valence band into the empty conduction band and/or the exciton levels, and observe 
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the spectrum of the subsequently emitted light coming from the eventual recombination of 
the excited electrons and the holes created in the valence band. The photoluminescence 
spectrum is proportional to the exciton/EHP optical oscillator strength, which, at first sight 
(i.e. without incorporating the Sommerfeld factor effect associated with the electron-hole 
Coulomb interaction), is expected to have an cu -1 / 2 singularity at the band gap energy due 

to the E~ x l 2 divergence of the ID electron density of states |L8| at band edge. However, 
this ID plasma band edge singularity is known to be strongly suppressed by the excitonic 
Coulomb correlation effect JXOj so that the main peak observed in the experimental pho- 
toluminescence spectra should result from the excitonic effect rather than the band-edge 
singularity of the noninteracting electron-hole plasma. The most striking experimental ob- 
servation in the recent experimental studies of photoexcited QWR systems has been 
the finding that the exciton peak seems to be at an almost constant energy indepen- 
dent of the carrier density, i.e. independent of the laser pumping power. Thus the exciton 
peak seems to remain well-defined (and unshifted in energy) all the way from very low to 
very high photoexcitation density (~ 3 x 10 6 cm -1 jTJ) without any distinct signature of 
the expected insulator(exciton)-to-metal(EHP) Mott transition and the associated optical 
gain. The constancy of the exciton energy could, in principle, arise from an almost exact 
cancelation between the exchange-correlation induced shrinking of the nominal band gap, 
the so-called band gap renormalization (BGR), and the reduction of the exciton binding 
energy (with respect to the bottom of the renormalized band edge) due to the screening 
induced softening of the Coulomb interaction 0. Such an accidental cancelation between 
two distinct physical mechanisms (namely, BGR and screening suppression of exciton en- 
ergy) over a wide range of photoexcitation density needs to be theoretically established in 
a compelling way |J. In addition, combining this accidental cancelation explanation with 
the experimental fact of very a high Mott density (not yet seen experimentally) one may 
conclude that the BGR of ID electron-hole (e-h) system should be very weak in the high 
density situation, which is not consistent with the theoretical calculations upto now P,P^B],^B 



In particular, one must understand why there is no characteristic signature of the EhP m 
the luminescence spectra even at very high photoexcitation densities. One must be able to 
answer the question: where has the Mott transition gone ? On this issue, an important 
and unresolved problem for the photoluminescence experiment is that there is no reliable 
and direct way of estimating the photoexcited electron-hole density in such highly pumped 
QWR systems. The theoretical basis of the density estimation methods in the literature 
0, such as from the lineshape analysis of the spectrum, is usually not self-consistent and 
not appropriate in such high density strongly laser-pumped systems. Although we feel that 
the precise carrier density of the photoexcited QWR systems may not be known accurately, 
this issue does not pose any fundamental problem for our theory where the EHP density 
n = n e = nh is an input parameter. The problem arises only in trying a direct quantitative 
comparison with experiments p|-p[. 

From the theoretical point of view, the full many-body calculation in a high electron-hole 
density semiconductor system is complicated and has not been attempted before except for 
our own short letter published last year |J. The exciton mode is a solution of the Bethe- 
Salpeter equation (BSE) for the interaction vertex which, in the many particle situation of 
interest to us, should include self-energy and dynamical screening correlations. A complete 
or exact solution of the BSE is only possible in the dilute exciton limit when it reduces to 
a simple hydrogenic electron-hole bound state Schrodinger equation. Our interest in this 
paper is in the many-particle "exciton" state in the photoexcited semiconductor QWR sys- 
tem where self-energy correlations of simple electron or hole states and dynamical screening 
of the electron- hole Coulomb interaction vertex are important. A model of an electron gas 
with a single hole in a wire |17]] is not appropriate in our problem because a bound state 
always exists in any attractive potential in ID systems, which will trivially provide an over- 
estimate of the Mott density. We emphasize that both the quasi-particle self-energy and 
the dynamical screening of the electron-hole interaction vertex should be included properly 
(i.e. consistent with each other in a conserving approximation) in the BSE to obtain the 
correct description of the Mott transition. With the exception of our own short earlier re- 
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port || most other theoretical calculations use the static (Hartree-Fock) approximation or 
the quasi-static approximation |Tj| to the self-energy and a statically screened interaction 
vertex to solve the many-particle BSE and obtain the optical absorption/gain spectra. In 
these simpler approximations, where dynamical screening effects are neglected in an uncon- 
trolled way, the dominant excitonic peak has a weak red-shift (a few meV decrease) with 
increasing density upto a Mott density, n c , above which the excitonic peak completely disap- 
pears and the spectrum shows a shallow (and weak) gain region very similar to the behavior 
of the noninteracting EHP ||7| |iO| , |l^|T7|1 . Including the many-body dynamical screening || in 
the Coulomb interaction, the exciton peak stays essentially constant in energy (for n < n c ) 
and exhibits a pronounced gain spectrum (for n > n c ) , stronger than the quasi-static results. 
But the predicted Mott densities in the above theories (n c ~ 8 x 10 4 - 8 x 10 5 cm -1 ) are all 
below the experimentally estimated value (n c > 3 x 10 6 cm -1 ) — in fact, it is not clear if 
experimentally the transition to the EHP has ever been observed even at the highest pho- 
toexcitation densities. It is in general hard to include the many-body effects appropriately 
in a calculationally tractable model over such a wide range of density (over at least 4 orders 
of magnitude in n), from the weak coupling dilute exciton gas system to the strong coupling 
EHP regime. 

In this paper, starting with the realistic Coulomb interaction in ID T-shaped QWR 
systems, we first evaluate the single particle self-energy for both electrons and holes in the 
dynamical plasmon pole approximation (PPA) within the so-called GW scheme (i.e. in the 
leading order dynamically screened interaction) to obtain the electron and hole renormal- 
ized Green's function. This self-energy calculation which by itself does not contain any 
direct excitonic effects, gives us the BGR or the reduction of the nominal band gap due to 
exchange-correlation. For comparison, we also calculate the BGR obtained by the quasi- 
static calculation in both static random phase approximation (RPA) and static PPA in 
this paper. We then derive analytically the effective electron-hole (e-h) interaction vertex, 
V e ff(k,u>), which includes consistently the electron-hole-plasmon coupling with the external 
photons within our dynamical GW-RPA-PPA approximation scheme. We use two different 
methods to study the excitonic properties: one is a variational approximation on an effective 
exciton Hamiltonian fL8 , which depends on the carrier density; the other technique is to 



solve the dynamical BSE by treating both self-energy renormalization and vertex correc- 
tion (arising from the Coulomb interaction) on an equal footing (within our plasmon pole 
approximation scheme), obtaining the optical absorption spectra. Both calculations are car- 
ried out over a wide range of e-h density from n = 10 2 cm -1 to n = 10 6 cm -1 at finite 
temperatures under the quasi-equilibrium condition, i.e. the e-h density is assumed to be a 
constant parameter for each density calculation (and n = n e = rih). While our dynamical 
BSE calculation includes exciton and EHP effects equivalently and is directly capable of 
providing the Mott density n c through the analysis of the absorption spectra, the varia- 
tional exciton energy has to be compared with the BGR calculation in order to ascertain 
the Mott transition — in particular, the merging of the effective variational exciton with 
the renormalized band edge is taken to be the signal for a Mott transition. We find that the 
absorption peak obtained from solving the dynamical BSE survives with very large broad- 
ening well above the critical density n c estimated from the variational approximation, and 
no optical gain (negative absorption) regime shows up in the spectra even at the highest e-h 
density. This implies the non-existence of Mott transition in ID electron-hole systems. This 
striking result may be physically understood as arising from the fact that the quasi-particle 
picture underlying the conventional Fermi liquid model fails in high density ID systems 
due to strong inelastic scattering by plasmons, associated with the generic non-Fermi liquid 
behavior in ID systems. In fact, in ID systems there is no conventional EHP because there 
are no single particle excitations in an interacting ID systems. This non-existence of single 
particle excitations or quasi-particles also leads to a breakdown of the conventional exciton 
picture — a quasi-electron and a quasi-hole bound pair. We will discuss this point in more 
details later in this paper. 

This paper is organized as follows: in Sec. II we present and discuss the theory we 
use in various parts of our calculations, e.g. the realistic Coulomb interaction in ID T- 
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shaped QWR system, the single particle self-energy calculation in the single loop PPA- 
GW approximation, the different approximation schemes used for screening the long-ranged 
Coulomb interaction, the dynamical Bethe-Salpeter equation approximations in our theory, 
and the effective exciton Hamiltonian used in the variational calculation. In Sec. Ill we 
show our results for the density dependent exciton energy in the variational method and for 
the excitonic optical properties from the solution of BSE. In Sec. IV we conclude with a 
discussion and a summary of our results. 



II. THEORY 

We use the two-band (one conduction band and one valence band) model to study the 
ID electron-hole system, neglecting higher subbands and the degenerate valence bands. We 
also consider the photoexcited quasi-equilibrium regime where the e-h density is assumed 
to be a constant (in time) so that the Hamiltonian of such a ID electron-hole system can 
be expressed as (in the effective mass approximation and assuming purely parabolic band 
dispersion; we take h = 1 throughout): 

+ i"E [VcMtii-qcl+qCk'Ck + V cM {q)d\_ q d\ l+q d kl d k - 2V c>eh {q)c\_ q c k d\, +q d k ?\ , (1) 



k,k',q 



where c k (c\) and d k (d\) are the annihilation(creation) operators for conduction band elec- 
trons and valence band holes respectively (we will not explicitly show the spin index in 
summations throughout this paper although spin is included in our calculations), and m e /h 
are the electron/hole effective masses. E® is the direct band gap between the top of the 
valence band and the bottom of the conduction band, taken to be 1550 meV for the 
GaAs/Ali-zGa^As QWR system in all our calculations. There are three different Coulomb 
interactions entering the Hamiltonian: electron-electron (V^ ee (g)), hole-hole (V c ^h(q)), and 
electron-hole (K,e/i(<?)) interactions. The first two give rise to the electron and hole quasi- 
particle self-energies and the other one, the electron-hole interaction, produces the exciton 
bound state. One should note that if we neglect the self-energy correction and also dy- 
namical screening effect (i.e. the low density limit of a Wannier exciton), the Hamiltonian 



of Eq. (1) leads to a ID hydrogen atom problem for the Wannier exciton, which in 
ID always has a bound excitonic state even for an arbitrarily weak electron-hole (attrac- 
tive) interaction. Using a model of an electron gas with a single hole therefore gives rises a 
very high Mott density estimate (even if V c ^ e h is statically screened), which is a reflection of 
this ID bound state property. We address both the many-body self-energy effect and the 
electron-hole excitonic binding effect on an equal footing in the theory, which we accomplish 
by using the dynamical Bethe-Salpeter equation as described below. 



A. Coulomb interaction in QWR 



The realistic (bare) Coulomb interaction in ID QWR is obtained by taking the expecta- 
tion value of the 3D Coulomb interaction over the electron wave function along the transverse 
directions (y and z axes, see the inset of Fig. 1) of the wire. After Fourier transformation 

I3JI5IJ 



along the 
element: 



ID wire direction (x), we have 



for the Coulomb interaction matrix 
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^2 yoo POO /■! 

= — / dy dy f dz dz' / 

£q J—oq J -co J — 



00 ^x 2 + (y - y') 2 + (z — z') 2 



dydy' dz dz'\Uy, z)\ Z \Hy'> z')\ 2 K (q^/(y - y'Y + (z - z'Y), (2) 



2e 2 f°° 

Eq J — 00 

where <pi(y,z) is the QWR confinement wavefunction for the lowest eigenstate of electrons 
(i = e) or holes (i = h). Their exact forms depend on the geometry and the detailed nature 
of confinement for the QWR system. Kq(x) is the zeroth-order modified Bessel function of 
the second kind |20| which diverges logarithmically when x goes to zero (i.e. in the long 
wavelength limit). 

Following the experimental system of Ref. j2], we use T-shaped QWR parameters to 
numerically calculate the ID Coulomb interaction via Eq. (2). To simplify calculations (and 
also to have some analytical control) we use the following two approximations in evaluating 
the wavefunction <pi(y,z): (i) we assume the confinement potential for both electrons and 
holes to be infinitely deep, i.e. both electrons and holes are completely confined by the 
2D T-shaped potential well, so that the wavefunctions of electrons and holes are of the 
same form, independent of their effective mass difference. Consequently the three different 
interactions (K jee , V c ,hhi an d V C)e h) become the same, denoted by V c throughout this paper. 
This simplifying approximation is justified by the detailed work of Ref. ||, as mentioned 
in the Introduction, (ii) Instead of numerically solving the complicated 2D Schrodinger 
equation to get the ground state single particle wavefunction |j5,12,T3| (which is not the 



focus of our interest), we simply approximate (p(y, z) to be the product of two single- variable 
functions, £(y) and ip(z) (i.e. <fi(y,z) ~ £(y)ip(z)), and assume plj a simple and reasonable 



approximate model form for £(?/) and ip(z) through the following exponential formulae: 

W z 

where W z and W y are the full-plane (x — y plane) QW width and the half-plane (x — z 
plane) QW width respectively. Eqs. (3) and (4) have the maximum electron/hole density 
at y — 0, z — W z /2, with three branches of exponentially decaying density along ±y and 
+z directions (see the inset of Fig. 1). The exponential decaying lengths or confinement 
sizes are W y /2 and W z /2 in y and z directions respectively, and thus in our model of the 
T-shaped QWR the effect of wire geometry on the Coulomb interaction is entirely contained 
in the effective "wire size" W y and W z , which are the confinement parameters of our model. 
This approximation greatly simplifies the calculation of the realistic Coulomb interaction in 
Eq. (2) and makes our BSE calculations tractable. We believe that our QWR confinement 



model, as defined in Eqs. (3) and (4), to be quite reasonable 21]. For example, the exciton 
binding energy calculated in this approximation is 18.2 meV tor W y = W z = 7 nm wire, 

very close to the quoted experimental value, 17 meV, for the same wire size 0. The small 
over-estimate (about 7%) is expected because of the assumption of infinite confinement 

energy and the strong e~ y localization of £(y). In more accurate numerical treatments the 
confinement is weaker than in our model, leading to a lower binding energy in the QWR 
system. In Fig. 1 we show the calculated V c (q) from Eqs. (2)- (4) for different wire sizes. 
We assume only one (the ground) electron and hole subband in the conduction and valence 
band respectively. 
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B. Absorption spectra 

In order to study the excitonic effect on optical properties of ID photoexcited electron- 
hole systems in semiconductor QWR structures, we calculate the dynamical (photon fre- 
quency dependent) absorption coefficient, a(u), and refractive index, n{ui), which are related 
to the long wavelength dielectric function, e(q — > 0,u>), by the following formula 

B ( w )+i£^= e ( w )V», (5) 

where c is the vacuum light velocity. The dynamical refractive index, n(u), is therefore given 
in terms of e(u) by, 



n(u) = J- \Ree(iu) + (Ree(uj) 2 + lme(uj) 2 ) 1/2 } (6) 



, 2 

and the absorption coefficient, a(w), is given by 



«H = ( \ ■ 7 
n{u))c 



Using the linear response theory flqJZl , the dielectric function of the ID e-h system is 
expressed as 

Aire 2 

e(w) ~ Eoo - — — r vc (k)r* c (k')G g ^o{k, k', u), (8) 
AL k,k' 

where the retarded pair Green's function, G q (k, k',u), is 

G q (k,k',u) = jr° e <wt ([d_ fc (t) Cjfe+fl (*) > 4^(0)^(0)])^^ (9) 

and A = W y W z is the cross sectional area of the QWR. In these equations q is the center of 
mass momentum of the exciton which is set to zero (and hence not shown explicitly) in all 
our calculations below. r vc (k) is the dipole matrix element, which can be simplified in the 
effective mass approximation ||18j| : 

where the reduced mass m = m e mh/(m e + m^) and 



M(k) 



( 1 + ^ei) 1 - 



By introducing a new function 

Q{k,u) =J2M{k')G(k,k',uj), (12) 
k' 

the dielectric function in Eq. (8) can be expressed to be 

eH ^ ^oo - E M(k)Q(k, W ), (13) 

g k 

and the dynamical function Q(k,u), which is essentially a two-particle Green's function, 
satisfies the Bethe-Salpeter equation described below. 
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C. Bethe-Salpeter equations 



For the results to be presented in this paper, the many-body exciton is given by the so- 
called Bethe-Salpeter equation |18| for the two-particle Green's function shown diagrammat- 
ically in Fig. 2(a), which corresponds to a rather complex set of two-component (electrons 
and holes) coupled non-linear integral equations which must be solved self-consistently with 
the bare interaction being the Coulomb interaction in the QWR geometry. These equations 
are notoriously difficult to solve without making drastic approximations, and in fact have 
never before been solved in the literature in any dimensions (except for our own short re- 
port earlier H]). In carrying out the full many-body dynamical calculations for BSE we are 
forced to maSe some approximations. Our most sophisticated approximation uses the fully 
frequency dependent dynamically screened electron-hole Coulomb interaction in the single 
plasmon-pole approximation, which has been shown to be an excellent approximation |23| 



to the full random phase approximation (RPA, see Fig. 2(c)) for ID QWR system. For 
the self-energy correction we use the single- loop GW diagram shown in Fig. 2(b) with the 
screened interaction given by PPA. Ward identities then fix the vertex correction, entering 
Fig. 2(a), to be the appropriate ladder integral equation. 

For convenience, we use the finite temperature imaginary time Matsubara frequency 
Green's function formalism in our analysis. The bare electron-hole two-particle Green's 
function without any e-h interaction is 

G°(k : k', z, fi) = G e (k, Q - z)G h (-k, z)5 k) y, (14) 

and it corresponds to the two separate Green's function lines of electron and hole in Fig. 
2(a). For each particle line, we have 

-«*-*<*«>+«• (! = e '' l) (15) 

where e e ^ = k 2 /2m e + Eg and e^j. = k 2 /2rrif l are the bare (noninteracting) band energies 

for electrons in the conduction band and for holes in the valence band respectively; [a is 
the chemical potential and z) is the self-energy (for a complex frequency, z), which 
we will calculate later within GW approximation. In order to avoid the multi-pole (and 
any possible branch cut) structure in Gi(k, z), we approximate z) by the momentum- 
dependent band gap renormalization, Aj(fc), which is related to the self-energy through the 
self-consistent Dyson's equation: Aj(ft) = Sj(fc, e^k + Aj(fc) — //j), i.e. Ai(k) is the so-called 
quasi-particle on-shell self-energy. However, Aj(&) can be well approximated |24| by trun- 
cating this equation at the first nontrivial order, i.e. Ai(k) = Ej(&,£gfc — which should 
be reasonably valid in our calculations below. Therefore we have the following electron/hole 
single pole Green's function, 

Gfa z) 1 , (16) 

z - e itk - Aj(fc) + ^ 

for later calculations in this paper. The details of calculating the self-energy Sj(fc, z) within 
the GW approximation are discussed in the Sec. II-D below. 

The Bethe-Salter equation in Fig. 2(a) could be read as (with (3 = l/k B T, where T is 
the temperature) 

G(k, k', z, Q) = G°(k, k',z,n)x I 1 + i v s( k - k "> z - z')G(k", k', z', Q) ] . (17) 

\ P k",z 

Putting Eqs. (14)-(16) into Eq. (17) we get 



S 



(ft - £ e , fc - efc,_ fc - A e (A;) - A fc (-A;) + /i e + fih)G(k, k', z, ft) 

= (G e (fc, ft - z) + G h (-k, z)) 5 k>k > x M + 1 ^ V s (k - fc", 2 - z')G(k", k', z', ft) ] . (18) 

\ P k",z' 

This equation, however, is not of closed form and is difficult to evaluate since it is a rather 
complex multidimensional singular integral equation. We therefore have to use an addi- 



tional simplifying approximation first introduced by Shindo |T8|J25|J26|| , where the two parti- 
cle Green's function, G(k, k', z, ft), is replaced by a simple pair Green's function G(k, k', ft) 
(whose retarded function yields via Eq. (9) directly the optical dielectric function): 



G(k, k', z, ft) ~ 

G e (k,Sl-z)+G h (-k,z) 
-^E z (G e (k } n-z)+G h (-k,z)) 



G(k,k',n), (19) 



where 



G(M',ft) = -^E G (MUft), (20) 

and 

-^(G e (k,Q-z)+G h (-k,z)) = 

1 - n e (^ k ) - n h (£ h - k ). (21) 

Here ^ = e^ k + Aj(/c) and ni(^ k ) is the fermion momentum distribution function, 
( e P(Tte£i,k-m) _|_ which keeps the electron and hole density constant by adjusting the 

chemical potential, fii, to satisfy the correct density constraint, / \dk / V) rii(£i t k) = n. Note 
that the approximation defined by Eq. (19) follows from the exact BSE in a statically 
screened Coulomb interaction i.e. if the frequency dependence of the effective dy- 
namically screened interaction is neglected. We expect the Shindo approximation to be a 
reasonable approximation in our dynamical calculation below, because the dynamical screen- 
ing effect contributes mostly to the correlation energy, whose real part is dominated by the 
(static) Hartree-Fock exchange energy in the high density region |19| (while the imaginary 
part of the correlation energy plays an important role in our calculations below). We have 
not been able to find a tractable way of solving the dynamical BSE without making the 
Shindo approximation. 

Using Eqs. (15), (19)-(21) in Eq. (18), we then have the following effective Bethe- 
Salpeter equation for the pair Green's function G(k, k', u) (after the analytical continuation 
ft — > to + id — /i e — jJh)'- 

G(k,k',u) =G°{k, k',u) 

x , (22) 

where G° and the dynamically screened effective electron-hole interaction, V e ff, are expressed 

G°(k, k', co) = \- n ^e,k)-n h ^ h k ) 

uo + ib - e e , k - e h k - A(k, to) 



as 
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and 



V eff (k, k',u) 



I 

J, 



E 



G e (k, n - z) + G h {-k,z) 



1 - ™e(£e,fc) - n h (£, h - k ) 

G e (k',n-z') + G h (-k',z'y 



V s (k-k',z- z') 



(24) 



£l=UJ —fle—fjtfr -\-iS 



1 - n e (£,e,k>) - n h (£, h - k >) 
The effective BGR, A(k,u), is given by 

A(k,tj) = [(1 - ™e(£ e ,fc'+<z) - n h (£ h ,-k>)) x k',q,uj) - V C (A; - fc')] <5 SS ' 

fc' 

= - EWW + «a(&,-*')]K//(*, + ?» ^) - V d k - (25) 

k' k' 

In Eq. (25) the self-energy term, (n e + rih)V e ff, and the vertex correction, V e ff — V c , are 
treated on an equal footing. G°(k, k' ,u) in Eq. (23) is the electron-hole pair Green's 
function with self-energy correction but without electron-hole attractive interaction, which 
is now replaced by the dynamically screened effective interaction, V e ff(k, k',u), in the BSE, 
Eq. (22). If we neglect dynamical effects in V s (k,z) (as in the static or the quasi-static 
approximation described below), then V e ff(k,k',uj) = V s (k) according to Eqs. (24) and 
(21). In the following section, we will discuss the use of different screening models to evaluate 
V e ff(k,k',u) (and BGR, A e /h(k), through the screened GW approximation) in calculating 
the absorption spectrum by solving the BSE. 

Combining the Bethe-Salpeter equation, Eq. (22), for G(k, k',uj) with Eq. (12), we have 
the following equation for Q(k,u): 



Q(k,u) = Q (k,uj) x 1 



M{k) 



Y,V eff (k,k',uj)Q(k\ 



(26) 



for Q (k,uj) = J2k' M(k')G°(k, k', u). Once Q(k,u) is obtained by solving the integral 
equation, Eq. (26), which is also a BSE, it is straightforward to calculate the absorption 
and gain spectra from the dielectric function, e(u), through Eq. (13). 



D. Self-energy, BGR and screening in QWR 

In order to solve Eqs. (22) to (26) for the Bethe-Salpeter equation, we have to use a 
screened interaction, V s (k,z), in Eq. (24) to get V e ff and also to get the single particle 
self-energy, Sj(A;, z), in the Green's function of Eq. (16). In this section, we discuss and 
compare both the quasi-static approximation and the dynamical (PPA) approximation in 
the screening calculation. For convenience, we first discuss the self-energy part and then the 
screening effect. 

In the GW approximation which is the leading-order self-energy in the screened inter- 
action expansion, the self-energy is calculated in the single-loop diagram composed of a 
noninteracting particle line and a screened interaction line (Fig. 2(b)). Using static screen- 
ing in the interaction line (i.e. V s (k,z) = V s (k,0)), we get a screened exchange self-energy 
term only, and all higher order screening effects to the correlation energy are neglected. 
This approximation (named static approximation) is therefore too simplistic to give correct 
results H, although it has been extensively employed in excitonic calculations because of 
its simplicity. An improvement to the static approximation is the quasi-static approxima- 
tion [jy| , which neglects the recoil energy during the scattering process so that no dynamical 
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frequency inside the screened interaction potential shows up. This approximation produces 
an extra constant Coulomb-hole term (the second term of Eq. (27)) in the self-energy in 
addition to the screened exchange self-energy of the static approximation, so that the full 
expression for the BGR in this quasi-static approximation becomes 



Ai(k) 



E 

k' 



-V s (k - k/)m{e iik ) + -(V s (k') - V c {k')) 



(27) 



where V s (k) = V s (k,u = 0) = V c (k)/e(k, uj = 0) is the statically screened Coulomb inter- 
action, which could be analytically derived either from RPA (using Eq. (29) below) or 
PPA (using Eq. (30) below) ||23|| . In our paper, the former is called quasi-static-RPA and 
the latter named quasi-static-PPA. Note that Aj(fc) in Eq. (27) is pure real, i.e. without any 
imaginary part of the self-energy or inelastic broadening effect, so that the quasi-particle as- 
sumption for the Landau-Fermi liquid is completely satisfied in this approximation with an 
infinite quasi-particle life time. It is well-known, however, that the quasi-particle assumption 
breaks down in ID (unlike in 2D or 3D) electronic systems, with a generic non-Fermi liquid 
behavior ||27|| . For the purpose of comparison we still use this approximation to calculate 
the ID optical properties in order to compare with the full dynamical calculation results 
and to study the quantitative validity of this widely-used quasi-static approximation both in 
the higher dimensional systems []18 , |j2| and in the ID system |7|JT0|] in the literature. In Fig. 
3(a) we show the conduction band energy, £° k — E ° = e e ^ + A e (k) — E®, in the quasi-static 
PPA for different electron densities. The band gap renormalization is almost a wavevector 
independent rigid shift in the quasi-static approximation. 

For the self-energy Sj(fc,u;) calculated in the one-loop GW approximation with dynam- 
ically screened interaction, we have 



z) = -- J2 v s(k - k', z - z'^k', z') 

P k',z> 

v ^- k ') C (k> z') 



(28) 



where we can use either RPA or PPA (which is an excellent approximation to RPA |2^j) to 
calculate the dynamical dielectric function, s(k, uj). For zero temperature RPA, e(k, uj) is 
obtained by including the non-interacting polarizabilities of electrons (n^fc, uj)) and holes 
(U° h (k,uj)) H: 



e{k,uj) = 1 - V c (k)U° e (k,uj) - V c (k)U° h (k,uj) 

mi , \u 2 — [(k 2 /2mi) — kv Fji } 2 
- [(fc 2 /2mi) + kv F ^ 2 



i=e,h nK 



(29) 



where VF, e /h is the (Fermi) velocity of electrons/holes at Fermi momentum in the conduc- 
tion/valence band. In this paper we will use RPA in only calculating the quasi-static screen- 
ing via Eq. (27) by setting uj = in Eq. (29), not in the full dynamical BSE (Eq. (26)), 
because the pole structure (and branch cut properties) of screened interaction, V c (k)/e(k, uj), 
in the full dynamical RPA is too complicated to deal with in the frequency summation of 
Eq. (24) and in the integral equation of Eq. (26). In the dynamical PPA, however, the 
dielectric function e(k, uj) is defined by the following expression where screening is induced 



by a single (plasmon) pole satisfying the corresponding /—sum rule [23 



e(k, uj) 



1 + 



OJ 



pi 



(A) 



[uj + isy 



OJ 



2 ' 



(30) 
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where uj p i(q) = JnV c (q)q 2 /m is the ID plasmon oscillator strength and LO q is the effective 



plasmon frequency given by a simple formula fl4|J18 



UJ 



pi 



9 + — + 

mK 



Am 2 



(31) 



where k is the inverse screening length. It has been shown that PPA is a very good ap- 
proximation to RPA in ID systems, where plasmon excitations dominate the single particle 
excitations [15|.|23|]. The great virtue of the single-pole PPA for our theory is that it makes 
our calculation of V e ff in Eq. (24) tractable because the integral equation in frequency be- 
comes simple. In the PPA the self-energy of electron (i = e) or hole (i = h) can be expressed 
as a sum of the usual exchange or Hartree-Fock energy, T,f x (k), and the correlation energy, 

k> 

, x _ V c (k')uj 2 (k') 
Zr(k,u;) = Y: o P x 



n B (uj k >) + rii{e i)k+kl ) n B (uj k >) + 1 - n^e^+y) 

UJ + UJy - Si^k'+k -ij UJ—UJy — Ei^k'+k + h 



(32) 



for the plasmons; 



where n B {uJk) is the bosonic momentum distribution function, (er Uk — 1 
7 is a small phenomenological damping term incorporating impurity scattering and all other 
possible broadening process (see the discussion in Sec. III-B). From Eq. (32) we see that the 
dynamical effect as well as the imaginary part of Ej(fc, uj) arises entirely from the correlation 
energy (and is absent in the static (Hartree-Fock) or quasi-static theory). This will play an 
important role (which is crucial in ID) in our following calculations. Figs. 3(b) and 3(c) 
show the real and imaginary parts of electron energy, £ efc — E° g = e e fc + S e (/c, £ e>k — /i e ) — E®, 

taking into account the dynamical PPA self-energy renormalization. 

Defining the on-shell self-energy to be Aj(/c) = — /i e ) where i = e,h, the 



imaginary part of A e (k) is proportional to the electron inelastic-scattering rate [15] arising 
from electron-electron interaction, which is very small when k is below some threshold 
momentum k c . For k > k c a new collective mode scattering channel opens up in which 
electrons lose energy by emitting plasmons. At zero temperature and for zero impurity 
scattering (clean system limit), it can be shown that the inelastic scattering rate diverges 

as (k — k c )~ 1 / 2 when k approaches k c from above in ID [15]. Note that this divergence in 
ImA e /ft(A;) also exists in the RPA calculation |0|, and is therefore a characteristic of the 
interacting ID system in the dynamical GW approximation, causing a gap to open up at 
k = k c in the real part of the self-energy as shown in Fig. 3(b). The existence of this gap in 
BGR (or the divergence in ImA e // l (/c)) reflects the breakdown of the quasi-particle picture 
in the ID electron system [27j within the perturbative GW approximation. An interacting 
ID electron system is known to be better described by the Luttinger liquid (LL) model 
than the Fermi liquid model due to the strong plasmon scattering effect arising from the 
limited phase space in ID. A Luttinger liquid, in contrast to a Fermi liquid, does not have 
any discontinuity in its momentum distribution function, and does not, therefore, have any 
true quasi-particles. The existence of a Luttinger liquid is a purely nonperturbative effect 
of interaction and happens in ID even for arbitrarily weak electron-electron interaction. We 
therefore cannot get a true Luttinger liquid within our perturbative GW approximation, 
but the opening of the gap in the real part of the self-energy (or equivalently the divergence 
in the imaginary part of the self-energy) is the perturbative signature of the breakdown 
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of the Fermi liquid picture. At finite temperature and for finite impurity scattering, the 
single particle properties calculated in the ID Fermi liquid model via the dynamical GW 
approximation are similar to the results calculated in the Luttinger liquid theory Therefore 
we believe that the strong inelastic scattering shown in Fig. 3(c) qualitatively reflects the 
LL character of ID systems, and our self-energy calculation is qualitatively correct for our 
purpose of calculating excitonic optical properties. Our inclusion of the strong inelastic 
scattering by plasmons catches some essential aspects of ID phase space restriction, which 
eventually leads to the nonperturbative formation of a ID Luttinger liquid, which is beyond 
the scope of this work. 

We evaluate the effective interaction, V e ff in Eq. (24), by using the same PPA approxi- 
mation and obtain: 



V ef f(k, k', u) — V c (k - k') 



1 + 



N eh (k) N eh (k 



—Xeh(k,k',uj) 



(33) 



where N eh (k) = 1 — n e (^ e fc ) — rih(£h,k) an d Xeh is given by the the following complicated 
formulae containing eight different terms associated with various dynamical processes in the 
ID e-h system. 



Xeh(k, k',u) 



u 2 pl {k - k') 



2cu 



x 



k-k 1 



-(1 + n B {uk-ki))n e (ie,k) + n B {uJk-k')ne{£e,k>) + w e (g e ,fc)w e (g e , fc /) 
-n B {oJk-k')n e {ie,k) + (1 + n B {uJk-k'))n e {ie,k') ~ n e (Ce,k)n e (Ce,k') 



ie,k — £e,fc' + Wk-k' 

n B (uk-k'))nh(€h,-k) + n B (uj k - k >)n h (i h _ fc /) + n h (£ h - k )n h (£ h _ fc < 



ih-k — ih,-k' — k'fc-fc' 

-n B (u k - k ,)n h (£ h - k ) + (1 + n B (vk-k'))nh(th,-k') - ,-k)n h (Ch -k') 



ih-k — ih-k' + ^>k-k' 

n e {ie,k)n h {ih,-k') + (1 + n B (u k _ k ,)){l - n e (^ k ) - n h (£ h _*,)) 



-n e (^e,k)n h (Ch-k') + n B (u k - k >)(\ - n e (£ e<k ) - ra/, (&_#)) 



UJ + 27 + /I e + fi h - £ eifc - £ h + UJ k _ k , 

n e {ie,k')n h {ih-k) + (1 + n B {u k _ k ,))(l - n e (^ k <) - n h (£ h - k )) 



-ne(€e,k>)n h (£h-k) + n B (uj k - k >)(l - n e (£ e>k ,) - n h (€h,-k)) 

UJ + 27 + fi e + fi h - £ e fe / - & _ fc + 



(34) 



where we use the same phenomenological damping parameter, 7, to broaden the resonant 
threshold energies in the denominators. The first two terms in the bracket of Eq. (34) 
describe the coupling of electron excitations with the plasmon, having the corresponding 
particle filling factors in the numerator and the resonance energy in the denominator. The 
third and the fourth terms describe the same plasmon coupling process for the holes. The 
first four terms are static and tu-independent in our approximation. The last four terms are 
dynamical and depend explicitly on uj. These last four dynamical terms describe processes 
which couple both electron and hole systems with the plasmon modes, and are extremely 
important in the dynamics of the photoexcited system. 
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E. Effective Hamiltonian and variational method 



Before solving the full dynamical Bethe-Salpeter equation, it is instructive to study 
the excitonic and the EHP effects separ ately by treating the influence of the EHP on the 
excitonic states as a perturbation ||18|J26|| . Using an effective Hamiltonian derived from the 
Bethe-Salpeter equation, we can variationally obtain the exciton ground state energy by 
minimizing the energy expectation value through an Is exciton trial wavefunction. The 
effective Hamiltonian treats the EHP as a perturbative effect and is written as H vp i(uj n ) = 



H pp> + Hpp'i^n), where 




Hl= K + £r )Spp>-V c (p- P ') (35) 



is the Hamiltonian for the single electron-hole pair with an unscreened Coulomb interaction 
(similar to a ID hydrogen atom) and the perturbation H' is 

H pp'( U n) = A(p, uj n )5 ppl + V c (p - p') 

-(1 - fe(€e,p) ~ fh(£h,-p))V e ff(p,p',U n ), (36) 

for the nth eigenstate of energy u n . Here we can explicitly see the physical meaning of 
A(p,u) and V e ff(p,p',u) analytically derived in Eqs. (24) and (25). We expect that the 
wave function of the exciton satisfies the corresponding Schrodinger's equation in the low 
density limit, where the screening effect is negligible. Thus this exciton effective Hamiltonian 
approach may be a reasonable approximation to calculate exciton energies and wavefunc- 
tions. 

For the exciton trial wavefunction, <p n {p)i in the momentum space, we use the two pa- 
rameter variational wavefunction first introduced by Nojima to express the ID exciton 
ground state as 



Mp) = V *w TWTi > (37) 

where A and a are two independent (positive) variational parameters in our calculation. 
K\(x) is the first order modified Bessel function of the second kind. This variational bound 
state wave function has the following form in the real space, 



exp — Jix/X) 2 + cr 2 
o (x) = J , (38) 

2a\K 1 (2a) 

where one can see that the variational parameter A represents the exciton radius and a 
smoothens or broadens the center of mass wavefunction at x — 0. We do not study the first 
excited state wave function, <fii(p), in this paper because it is not particularly relevant to 
the Mott transition process we are interested in, although the variational technique can be 
adapted to study excited excitonic states M. 



III. RESULTS 



We first show the variational results because conceptually this is the simplest approach 
since it is based on an effective single particle Hamiltonian. We obtain the BGR and the 
exciton binding energy by variational method in both the quasi-static approximation and the 
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dynamically screened GW approximations for various photoexcited carrier (e-h) densities. 
The crossover between the exciton energy and the BGR gives us an estimated Mott transition 
critical density, n c , where the exciton bound state ceases to exist, and an insulator-to-metal 
transition occurs. The idea here is that at n c the exciton merges with the e-h continuum 
and is no longer a stable bound state. Finally we carry out the full Bethe-Salpeter integral 
equation solution by a matrix inversion method and obtain the absorption spectra and 
refractive index in a large range of e-h density (from 10 2 to 10 6 cm^ 1 ) to compare with the 
variational effective Hamiltonian results. Details are described below. 



A. Effective Hamiltonian result 

In Fig. 4(a), we show the calculated density dependence of the exciton ground state 
energy variationally obtained from the effective Hamiltonian method and the BGR (A e (0) + 
A/j(0)) calculated in both the quasi-static approximation and the dynamically screened GW 
approximation as described in the section II-D. Both RPA and PPA are used in the quasi- 
static calculation (Eq. (27)) for comparison whereas the full dynamical calculations are done 
only in PPA. The intersection between the exciton energy (dashed lines) and the BGR (solid 
lines) indicates the Mott transition, where the exciton merges with the band continuum and 
the system has a phase transition from an insulating exciton gas to a conducting EHP. Note 
that the variational method introduced in Sec. II-E loses its accuracy near the Mott density 
(and becomes essentially meaningless for n > n c ), because the variational energy expectation 
value, E(X, a) = (4>o(X, cr)\H° + H'(E(\, a))\<f>o(X, a)), has a very flat minimum region in the 
A — a space around n « n c , i.e. the exciton wave function is highly broadened, so that its 
minimum energy is hard to determine in such perturbation-based variational method. In Fig. 
4(b) we show the variationally calculated trial exciton Is (ground state) wavefunction, (f>o(x), 
for different exciton densities. The exciton density dependences of variational parameters, 
A and a, are also shown in the inset of Fig. 4(b). The sharp divergences of A and a at 
n c ~ 2 x 10 5 cm -1 indicate the derealization of exciton ground state wave function, a signal 
of exciton-to-EHP Mott transition. In Fig. 4(a) we terminate the variationally calculated 
exciton line (dashed) at n = 2 x 10 5 cm -1 and use the dotted line to represent the peak 
position of the absorption spectra obtained from solving the full dynamically screened BSE 
(discussed below) to continue the exciton line to higher densities upto 6 x 10 5 cm -1 . 

We can make the following comments about the results shown in Fig. 4(a): (i) For 
density below 10 4 cm -1 the exciton energy has only a few meV density-dependent red-shift 
in the quasi-static- RPA/PPA approximations and almost no shift (less than 0.5 meV) in the 
dynamical screening approximation. This shows the almost complete cancelation between 
the exchange-correlation induced BGR (a density-dependent shift) and the blue-shift of the 
exciton energy (due to screening) over a wide range of density. On the other hand, using 
the static screening (i.e. exchange energy only) approximation in the same calculation does 
not lead to this cancelation ||, showing that the experimentally observed constancy of the 
exciton energy as a function of the photoexcited e-h density is a dynamical effect, which 
may not manifest itself in simpler approximations, (ii) For an e-h density higher than 10 4 
cm -1 , the exciton energy in the quasi-static- RPA has a rather large red-shift until it merges 
with the BGR line smoothly at n c ~ 6 x 10 4 cm -1 , indicating a rather low density Mott 
transition in this system. On the other hand, the exciton energies calculated in both the 
quasi-static-PPA and the full dynamical PPA are almost constant upto n = n c ~ 3 x 10 5 
cm -1 , where the band continuum meets the exciton energy, (iii) In the full dynamical results 
obtained by solving the dynamical BSE directly, the excitonic absorption peak (dotted line) 
seems to survive even for densities higher than the "critical density", n c , at which the 
calculated exciton energy crosses the band continuum BGR line. This shows that there 
must be reasonably strong hybridization between the exciton and the EHP in the dynamical 
BSE (note that the dotted line in Fig. 4(a) is not from the variational calculation, but is 
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obtained from the BSE solution), so that the effective BGR, including excitonic effects in 
the BSE, is actually less than the result we calculate from Eq. (32) by adding the electron 
and hole self-energies without incorporating exciton effects. This also demonstrates the 
limitation of the quasi-static approximation and confirms the necessity of the full dynamical 
BSE calculation in the high density ID e-h system. 



B. Dynamical Bethe-Salpeter equation result 

In Fig. 5, we show our calculated absorption and gain spectra by solving the full Bethe- 
Salpeter (integral) equation in the quasi-static and the full dynamical screening approxima- 
tions for W y = W z = 7 nm wire at a low temperature of T = 10 K. The integral equation for 
the two-particle Green's function, Eq. (26) (or equivalently Eq. (22)), is solved by the ma- 
trix inversion method with maximum momentum upto k max = (tt/2) x 10 8 cm -1 (= lOOfcp 
for n = 10 6 cm -1 ). The poles of the dynamical screened interaction, V e ff, in Eqs. (33) 
and (34), together with the logarithmic singularity of the ID Coulomb interaction in the 
long wavelength limit, produce a multi-singular kernel with multiple momentum-dependent 
singularities which have never been solved before in the literature (except for our earlier 



work M), because the usual singularity- removal method is ineffective here JT^.[22|]. In our 
calculations presented in this paper, we use a rather large matrix (1500 x 1500 in a Gaus- 
sian quadrature for \k\ < k max ) in the matrix inversion method in order to get good overall 
accuracy, i.e. the same calculations using even larger (2000 x 2000) matrix size (which is ex- 
tremely time-consuming and not shown here) do not show any significant difference (within 
10%) in the whole absorption (and refractive index) spectrum from the results we present 
in Figs. 5-7. The broadening 7 used in our calculation is a phenomenological parameter 
which simulates in a simple manner the effects of all possible scattering and broadening 
processes not explicitly included in our theory. These are, for example, impurity and defect 
scattering, inhomogeneities in the system (e.g. fluctuations in the wire width), broadening 
associated with optical excitation process itself, and phonon scattering. We mention that 
inelastic plasmon scattering life time effects are explicitly included in our theory. Note that 
7 should be small compared with the bare excitonic binding energy (~ 10 — 20 meV in GaAs 
semiconductor QWR systems), and as long as 7 is small, its precise value has no qualitative 
effect on our conclusions and results. We typically choose 7 = 0.5 meV in our calculations. 

In Fig. 5(c) and 5(d) we show the absorption spectra in the dynamical PPA for two 
different values of impurity scattering 7 (different by a factor of 2.5) to show that 7 does 
not affect the qualitative behavior of the spectra, but does control the linewidths of the 
absorption peaks as one would expect. Some important features of the optical spectra (cal- 
culated by solving the full BSE) shown in Fig. 5 are: (i) there are generally two absorption 
peaks in the low density (n < 10 4 cm -1 ) spectra of all the three approximations, one is the 
exciton ground state (Is) peak at about 1532 meV and the other one is the exciton first 
excited state (2s) at, for example, 1547.5 meV for n = 10 2 cm -1 (this peak is off the plot 
region in Fig. 5(a)). Note that this low density spectrum is almost the same in all the three 
different approximations, showing that the dynamical effect is not important in the low 
density region, (ii) When the density increases but is still less than 10 4 cm -1 , the exciton 
peak does not shift much (< 2 meV) with increasing carrier density in all approximations, 
indicating the constancy of the exciton energy, (iii) At higher densities (10 4 cm -1 < n < 10 5 
cm -1 ), however, the quasi-static- RPA result (Fig. 5(a)) shows some additional red-shift in 
the excitonic peak, consistent with the result shown in Fig. 4(a) which is obtained from the 
variational method. On the other hand, the excitonic peak positions in the quasi-static-PPA 
and in the full dynamical approximation are almost (density independent) constants in this 
region. A significant difference between the quasi-static-PPA and the dynamical calculation 
results, however, is that the exciton peak of the quasi-static-PPA results (Fig. 5(b)) has 
an almost constant oscillator strength, while the oscillator strength of the peak in the full 
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dynamical calculation results (Fig. 5(c)) decays at high density to about one-third of its low 
density value, (iv) For n > 10 5 cm -1 , both quasi-static- RPA and quasi-static-PPA results 
show negative absorption (gain) for the photon frequency below some critical value, u c , while 
the full dynamical result is still positive (i.e. no gain) with a weaker broadened peak upto 
n ^ 6 x 10 5 cm -1 or higher. In other words, we do not explicitly find the expected exciton 
(insulator) to plasma (metal) Mott transition when both the self-energy and the screened 
interaction are included dynamically in the full BSE theory upto a rather high e-h density. 
We believe that this behavior arises from the strong plasmon scattering effects in ID as 
discussed in Sec. IV of this paper. (Such strong inelastic scattering was not included in our 
earlier short report |J, leading to the appearance of a gain in the high density spectra above 
the Mott density.) In Fig. 6, we show the refractive index, n(u), calculated by solving the 
Bethe-Salpeter equation in both the quasi-static-PPA and the full dynamical approxima- 
tion for different photoexcitation densities. The calculated refractive indices in these two 
approximations are similar in structure. 

In Fig. 7, we show the calculated absorption/gain spectra obtained in both the quasi- 
static-PPA (Fig. 7(a)) and the full dynamical calculations (Fig. 7(b)) for the same wire 
width, W y = W z = 7 nm, but at a higher temperature (T = 100 K) for various densities. We 
find that the higher temperature (100 K) low density (n < 10 4 cm -1 ) absorption spectrum 
is almost the same as the corresponding lower temperature (T = 10 K) spectra in Figs. 5(b) 
and (c), while in the higher density region (n > 10 4 cm -1 ) the high temperature exciton 
absorption peak of the full dynamical calculation has a smaller red-shift in energy with a 
much larger broadening than the quasi-static-PPA result. The quasi-static-RPA result at 
such high temperature (not shown here) has an even larger red-shift and broadening. We 
mention that the gain in the absorption spectra of the quasi-static calculation at the lower 
temperature (Fig. 5(b)) is flattened and almost disappears in the higher temperature (100K) 
calculation results. 



IV. CONCLUSION 



In this paper, we theoretically study, using reasonably realistic Coulomb interaction, 
the excitonic optical properties of a ID QWR system by solving the many-body Bethe- 
Salpeter equation using a number of approximations, the most sophisticated one being a 
treatment of both the self-energy and the vertex function in the dynamically screened GW 
approximation. Our calculation is applied to the experimentally studied T-shaped GaAs- 
ALrGai-jrAs ID QWR systems for various densities and temperatures. We calculate the 
electron and hole self-energies in the one-loop GW approximation diagram using different 
screening approximations: the quasi-static-RPA, the quasi-static-PPA, and the dynamical 
(PPA) approximation. The quasi-static approximations give an almost rigid shift (the BGR 
effect) to the band energy (see Fig. 3(a)), and there is no imaginary part of the self-energy, 
i.e. the quasi-particle life time is infinite. This approximation may work well in 2D and 3D 
systems but fails completely in ID systems, because unlike in the higher dimensional systems, 
the electrons in ID system suffer very strong inelastic scattering effect by virtue of restricted 
phase space. Therefore the validity of the quasi-static approximation applie d to ID systems, 
which has been extensively used in many theoretical works |7,10Jl4,28,2^], is doubtful. In 
the dynamical calculation we find that the electron and hole oandgap renormalization has 
a gap opening up in its real part and a consequent divergent singularity in its imaginary 
part at k = k c (Figs. 3(b) and (c)), where the quasi-particle energy is transferred to the 
plasmon excitations due to very strong inelastic scattering by ID plasmons. Although this 
perturbative GW self-energy is "unphysical" due to the failure of the Fermi liquid model in 



the ID system [27], it still gives a rather good qualitative description of the single-particle and 



the collective mode properties (compared to the correct Luttinger liquid model), particularly 
at finite temperatures and for finite impurity scattering [15]. Our results in Figs. 3(b) and 
(c) reflect an important generic feature of ID systems: the quasi-particle excitation has a 
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very short life time (in fact, it does not exist) if the excitation momentum is higher than 
some value k c . This ID feature associated with Luttinger liquid properties of ID systems 
strongly affects the stability of ID excitons, the bound quasi-electron and quasi-hole pairs, 
as we can see from the calculated absorption and gain spectra (Fig. 5). 

In Fig. 5 we find that the quasi-static approximation, which excludes inelastic scattering, 
gives rise to a negative absorption (gain) region in the highly photoexcited system. The 
existence of gain means that the exciton state is saturated (fully occupied), and therefore 
manifests a spontaneous emission, rather than absorption. On the other hand, the overall 
positive absorption (no gain) spectrum found in the dynamical calculation (Fig. 5(c)) upto 

the highest density is caused by the large imaginary part of the electron/hole on-shell self- 
energy, lmA e / h (k) (see Fig. 3(c), and Eq. (32)), which is proportional to the inelastic- 
scattering rate and results from the energy scattering through plasmon channel. In other 
words, the excitons, composed of bound pairs of quasi-electrons and quasi-holes, are unstable 
due to strong inelastic scattering by ID plasmon excitations in the high density region. 
Consequently in the dynamical calculation, the exciton absorption peak is suppressed in 
strength and broadened in width as the photoexcitation density increases leading to stronger 
plasmon scattering. The absorption spectrum does not exhibit a negative (gain) region even 
in the high density regime because the quasi-particle EHP band continuum is so strongly 
inelastically scattered by plasmons that it is not a proper eigenstate (i.e. it decays) and 
is never saturated. The disappearance of the exciton line and the non-negativity in the 
absorption spectra (at the same time) in our dynamical calculation suggest that there is no 
insulator (exciton) to metal (EHP) Mott transition in ID systems, since both excitons and 
quasi-particles are strongly inelastically scattered by plasmons leading to neither of them 
being well-defined coherent states of the high density ID system. This result is consistent 
with the well-known non-Fermi-liquid properties of ID electronic systems, where the quasi- 
particle (and hence the exciton) picture fails. The quasi-static approximation, which ignores 
any plasmon effect and works well in 2D and 3D systems |2£| , does not work in ID systems 
because the ID excitation spectrum is completely dominated by plasmons. 

Another clue in support of the importance of plasmons in such high density ID e-h 
system comes from the temperature dependence of the absorption spectra shown in Fig. 7. 
Our results show that the high temperature (T = 100 K) absorption peak in the dynamical 
calculation (Fig. 7(b)) is suppressed and broadened so greatly that there is almost no 
spectral structure observed for n > 10 5 cm -1 , while the high temperature quasi-static- 
PPA result (Fig. 7(a)) still has a rather strong peak at the same density. This is because 
the plasmon excitation occupancy, whose energy distribution function, nsi^k), follows the 
Bose-Einstein statistics, depends strongly on temperature, leading to a qualitative difference 
between the T = 10 K and T = 100 K results in the dynamical calculation, while such 
plasmon dynamics is not included in the quasi-static calculation. This characteristic strong 
temperature dependence is also consistent with very recent experimental results ||. 

Based on our results and the discussion above, we propose that a crossover from a low 
density (essentially noninteracting) Fermi liquid to a high density interacting non-Fermi 
liquid is occurring in the optical spectra of the ID e-h system as the photoexcitation density 
increases (see Fig. 5(c)). In the low density limit, say n < 10 2 cm -1 , we have a dilute and 
noninteracting exciton system, whose absorption spectrum is independent of the many-body 
screening approximations we use — plasmons are just not that important in this regime. 
This shows that excitons in this situation are isolated quasi-electron and quasi-hole pairs, 
reflecting the validity of the quasi-particle picture in the effective noninteracting Fermi liquid 
model in the low density limit. In the higher density region, however, the plasmon effect 
on the quasi-particle self-energy becomes important, because the band curvature at k = kp 
is less for higher kp (i.e. higher density) and the relative importance of collective mode 
excitations (plasmons) is then strongly enhanced as in the Luttinger liquid model. Therefore 
the oscillator strength of the exciton absorption peak is then reduced and broadened in 
our dynamical calculation (Fig. 5(c)). When the density is roughly the nominal "Mott 
transition" critical density, n c , where the band continuum energy equals the exciton energy 
(see Fig. 4(a)), the plasmon excitation becomes so dominant that both exciton and band 
continuum (EHP) states become unstable, showing a crossover to effectively non-Fermi liquid 
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properties. We therefore do not expect to see the real Mott transition from an excitonic 
insulator to an EHP metal in ID highly photoexcited systems, in contrast to the results 
of previous theories. For an electron-hole plasma without any backward scattering in the 
usual Luttinger liquid model (no band curvature at all, which is unrealistic in our case), we 
can prove that gain in the optical spectra does exist below the Fermi energy at all densities 
with a complicated power-law singularity at the Fermi surface. Including the electron-hole 
attractive backward interaction (assuming a short-ranged interaction as in the so-called 
g— ology formalism), the electron- hole system most likely undergoes a charge density wave 
ground state transition with a mass gap in the elementary plasmon excitation ||30|| . While 
this scenario is consistent with our results, further work for the excitonic effect at energy far 
below the Fermi energy is still needed, because the regular Luttinger liquid model cannot 
include band curvature in an appropriate way in order to study the Mott transition at an 
energy level around the band edge. 

In reference to the experimental data, we note that our results from solving the dy- 
namically screened Bethe-Salpeter equation are in excellent qualitative and quantitative 
agreement with the recent experimental findings In particular, the effective constancy 
of the exciton peak as a function of the photoexcitedcarrier density as well as the possibility 
of excitonic absorption well into the high density regime (even for n > 6 x 10 5 cm -1 ) turn 
out to be characteristic features of the full dynamical theory (but not of the static and the 
quasi-static approximation). A full dynamical self-consistent theory as developed in this 
paper is thus needed for an understanding of the recent experimental results. Moreover we 
find that in our theory, the plasmon effect is crucial in the high density regime leading to 
the non-existence of any observable Mott transition in our calculation. This is consistent 
with recent experimental results [fflpj, which do not observe an actual Mott transition in 
the semiconductor QWR system even in the high photoexcitation density (~ 3 x 10 6 cm -1 ) 
regime. We emphasize that only our dynamical theory, and not the static or quasi-static 
approximation, is in agreement with the experimental results. 

In summary, our main accomplishments reported in this paper are the following: (i) the 
first fully dynamical theory of a photoexcited electron-hole system in semiconductors which 
treats self-energy, vertex corrections, and dynamical screening in a self-consistent scheme 
within a realistic Coulomb interaction-based Bethe-Salpeter theory; (ii) a reasonable quali- 
tative and quantitative agreement with the recent experimental observations of a constant 
(photoexcitation density-independent) excitonic absorption peak in energy, which in our 
dynamical theory arises from an approximate cancelation between the self-energy and the 
vertex corrections in the Bethe-Salpeter equation; (iii) inclusion of the plasmon effect in 
the quasi-particle self-energy calculation in our dynamical theory, leading to our theoretical 
proposal that no Mott transition should be observed in ID electron-hole systems (at least 
in optical experiments) even at very high photoexcitation density, i.e. there should be no 
optical gain region; (iv) instead, we suggest an experimentally observable crossover from a 
low density noninteracting Fermi liquid behavior (quasi-particle/exciton favored) to a high 
density interacting non- Fermi liquid behavior (no stable quasi-particles and excitons). A 
more precise and nonperturbative theoretical model for the high density ID electron-hole 
system is needed for future study — such a study should somehow incorporate both band 
curvature and Luttinger liquid behavior in analyzing the optical properties, although we 
believe that the qualitative features of such a theory are already contained in our work. 

This work has been supported by the US-ONR and the US-ARO. 
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FIGURES 
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FIG. 1. Theoretically calculated (from Eqs. (2)-(4)) ID Coulomb interaction in a T-shaped 
QWR system in momentum space. Results of different wire widths are calculated and shown 
together. In the inset we show the T-shaped intersection of two quantum wells in cross section. 




(c) 

FIG. 2. Many-body Feynman diagrams used in the paper with the single (double) solid line 
representing the bare (dressed) electron or hole Green's function and the single (double) wavy line 
representing the bare (dressed) Coulomb interaction: (a) the excitonic Bethe-Salpeter equation; (b) 
the single-loop self-energy (in the so-called GW approximation) defining the dressed Green's func- 
tion; (c) the RPA dressing of the Coulomb interaction (treated in the plasmon-pole approximation 
in our calculation). 
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FIG. 3. (a) Conduction band energy, — Eg, calculated in the GW approximation with 
screened interaction approximated by the quasi-static-PPA. (b) and (c) are respectively the real 
and imaginary parts of band energy calculated in the dynamically screened GW approximation for 
the same system as (a). The calculation is carried out in the symmetric T-shaped QWR system 
with W y = W z = 7 nm including finite temperature (T = 10 K) and finite (phenomenological) 
impurity scattering (7 = 0.5 meV) effects. 
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FIG. 4. (a) Separately variationally calculated exciton energy (dashed lines) and BGR of the 
EHP (solid lines) as a function of photoexcitation density, in three different approximations as 
indicated in the plot by different line widths. Note that when the density is larger than 2 x 10 5 
cm -1 , the variational method (introduced in Sec. II-E) fails to give a good exciton energy (see the 
text) and the dotted lines are the exciton peak positions of the corresponding absorption spectra 
by solving the BSE (Fig. 5). (b) The Is exciton ground state wavefunction obtained in the 
variational method through effective Hamiltonian (Eqs. (35) and (36)) in the dynamical (PPA) 
screening calculation for various electron- hole densities. Inset: the variational parameters, A and 
a, for the Is exciton ground state trial wavefunction with respect to the photoexcitation density 
in logarithm scale. When the density is near the Mott density (n c ~ 3 x 10 5 cm -1 ), both A and a 
increase sharply and the wave function becomes totally broadened. 
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FIG. 5. Calculated absorption and gain spectra for various photoexcitation densities by 
solving the Bethe-Salpeter equation in three different approximations for screening: (a) the 
quasi-static-RPA, (b) the quasi-static-PPA, (c) and (d) the full dynamical (PPA) calculation. The 
system parameters of (a)-(d) are the same as used in Fig. 3, except for the smaller 7 (=0.2 meV) 
used in (d). 




2 r 1 

1525 1530 1535 1540 1545 1550 
U (meV) 



26 



5 



2 



: (b) 


, r ,2 -1 - 

n= U cm 


io 4 : 




5.62X'0 4 : 




5.62X"0 5 : 




1Q 6 - 






7 





1525 1550 1555 1540 1545 1550 
CO (meV) 



FIG. 6. Calculated refractive index for various photoexcitation densities in both (a) the 
quasi-static-PPA and (b) the dynamical (PPA) approximation for interaction screening. 
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FIG. 7. Absorption and gain spectra obtained by solving the Bethe-Salpeter equation in 
(a) the quasi-static-PPA and (b) the full dynamical (PPA) approximation for screening at high 
temperature (T = 100 K). Other system parameters are the same as used in Fig. 3. 
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